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Abstract 

We consider solutions to the linear wave equation Ogcf) = 0 on a (max¬ 
imally extended) Schwarzschild spacetime. We assume only that the so¬ 
lution decays suitably at spatial infinity on a complete Cauchy hyper¬ 
surface E. (In particular, the support of (f) may contain the bifurcate 
event horizon.) It is shown that the energy flux of the solution through 
arbitrary achronal subsets of the black hole exterior region is bounded 
by C(uy^ -|- u^^), where v and u denote the infimum of the Eddington- 
Finkelstein advanced and retarded time of the subsets, and v+ denotes 
max{l,u}, etc. (This applies in particular to subsets of the event hori¬ 
zon or null infinity.) It is also shown that (fi satisfies the pointwise de¬ 
cay estimate 101 < in the entire exterior region, and the estimate 

|r0| < C^(l-|-|u |)“2 in the region {r > iZ}n J^(E), for any R > 2M. The 
estimates near the event horizon exploit an integral energy identity nor¬ 
malized to local observers. This estimate can be thought to quantify the 
celebrated red-shift effect. The results in particular give an independent 
proof of the classical result |0| < C of Kay and Wald without recourse to 
the discrete isometries of spacetime. 


1 Introduction 

The concept of a black hole is a central one in general relativity: Spacetime is 
said to contain a black hole when it admits a complete null infinity whose past 
has a regular future boundary. This boundary is called the event horizon and 
the black hole itself is defined to be its future; the past of null infinity is known 
as the black hole exterior. 

The simplest solutions of the Einstein vacuum equations of general relativity, 

= 0 ( 1 ) 
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containing black holes, the one-parameter Schwarzschild family of solutions, 
were written down in local coordinates m in 1916, but only correctly under¬ 
stood as describing spacetimes with black holes in the sense above, around 1960. 
Until that time, there were many arguments in the physics literature (e.g. cni) 
purporting to show that such solutions would be pathological and unstable. It 
was only when Kruskal demonstrated cni that the event horizon could be cov¬ 
ered by regular coordinates that its true geometric character became clear, and 
the problem of stability could be given a sensible and well-defined formulation. 

The Schwarzschild family turns out to be a sub-family of the two-parameter 
Kerr family which describe stationary rotating black holes. In its proper rigor¬ 
ous formulation, the problem of nonlinear stability of the Kerr family is one of 
the major open problems in general relativity.^ At the heuristic level, however, 
considerable progress has been made in the last 40 years towards an understand¬ 
ing of the issues involved. In particular, a very influential role was played by the 
work of R. Price m in 1972, who discovered a heuristic mechanism allowing for 
the decay of scalar field linear perturbations on the Schwarzschild exterior. The 
mechanism depends on the following fact, known as the red-shift effect, which 
had been understood previously in the context of geometric optics. Given two 
observers A and B, depicted below^. 



then if A emits a signal at a constant rate with respect to his own proper time, 
the frequency of the signal as received by B is infinitely shifted to the red as 
B’s proper time goes to infinity. 

For spherically symmetric solutions of the coupled Einstein-scalar field sys¬ 
tem, and more generally, the Einstein-Maxwell-scalar field system, the heuristic 
picture put forth by Price is now a theorem |^, and the red-shift effect de¬ 
scribed above plays a central role in the proof. Moreover, one of the results 
of |H], namely the decay rate 

101 + |a,0| < (2) 

along the event-horizon^, has important implications for the causal structure of 
the interior of the black hole: In [7|, it is shown using that, in the charged 

^In particular, it is conjectured that perturbations of Schwarzschild initial data should 
evolve into a spacetime with complete null infinity whose past “suitably” approaches a nearby 
Kerr exterior. 

^What is depicted is the Penrose diagram of a subset of Schwarzschild. See Sectionl2| The 
reader unfamiliar with these diagrams can refer to 1^ . 

^Here, u is a naturally defined Eddington-Finkelstein-like advanced time coordinate. 
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case, the spacetimes of |S] do not generically terminate in everywhere spacelike 
singularities as originally widely thought, but rather, their future boundary 
has a null component across which the spacetime metric can be continuously 
extended. In particular, this implies that the C°-inextendibility formulation ^ 
of strong cosmic censorship is false for the system considered. 

For the Einstein vacuum equations d) in the absence of symmetry assump¬ 
tions, results analogous to dCl seem out of reach at present. Clearly, a hrst 
step towards attacking geometric non-linear stability questions is proper under¬ 
standing of the linear theory in an appropriate geometric setting. This will be 
the subject of the present paper. Our main result is the following: 

Theorem 1.1. Let (f> be a sujficiently regular solution of the wave equation 

□g</. = 0 (3) 

on the (maximally extended) Schwarzschild spacetime {A4,g), decaying suitably 
at spatial infinity on an arbitrary complete asymptotically flat Cauchy surface 
E. Fix retarded and advanced Eddington-Finkelstein coordinates u and v on one 
of the exterior regions. For any achronal hypersurface S in the closure of this 
region, let F(S) denote the flux of the energy through S, where energy is here 
measured with respect to the timelike Killing vector field. Let = max{u, 1}, 
= max{u, 1}, and u+(iS) = max{inf 5 u, 1}, u+(5) = maxiinfg m, 1}. We 
have 

F{S)<C{{v+{S))-^ + {u+{S))-^). (4) 

(We also allows to be a subset of null infinity, interpreted in the obvious limiting 
sense.) In addition, we have the pointwise decay rates 

\(t)\ < Cvf^ in l^fiT+ycU+i^ 

|r(/)| < Cf^{l + \u\)-^ in {r > i? > 2M} n J+(E). (5) 

In the spherically symmetric case, the above result follows from a very special 
case of IH]. (See also m-) Decay for (f>, without however a rate, was first proven 
in the thesis of Twainy EH The uniform boundedness of (/> is a classical result of 
Kay and Wald m For the more general Kerr family, even uniform boundedness 
remains an open problem (see however EOl)- 

We should also note that, independently of us, a variant of the problem 
considered here is being studied by |3]. 

A statement of Theorem o in local coordinates will be given later. This, 
in particular, will explain the dependence on C, and the minimum regularity 
assumptions necessary. The reader wishing to penetrate deeper into this prob¬ 
lem, however, is strongly encouraged to learn the language necessary for the 
above geometric formulation. For neither the correct conditions on initial data, 
nor the desired statement of decay, have particularly natural formulations when 
stated with respect to Regge-Wheeler coordinates.^ The motivation for consid- 

^In particular, it is not correct to restrict to “compactly” supported initial data in Regge- 
Wlieeler coordinates, because r* = —oo corresponds to points in the actual spacetime. These 
issues are well known in the relativity community, and the reader should consult the nice 
discussion in m 
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ering decay as stated in 0 is that it is this formulation that has direct relevance 
both in the astrophysical regime, as well as for the fate of observers who enter 
the black hole region. For 0 applied to 5 = gives decay rates for the energy 
radiated to infinity (this is what is astrophysically observable), while applied to 
S — Ti.^, it gives decay rates for the energy thrown into the black hole (this is 
what concerns the observer entering the black hole). 

It is interesting to note that the decay rates 0^0 are sufficiently fast 
so as to suggest that the picture established in [7] may remain valid in the 
absense of symmetry assumptions, in particular, the existence of a marginally 
trapped tube which becomes achronal and terminates at and a weak null 
singular boundary component® to spacetime, emanating from i’*', across which 
the spacetime can still be continuously extended. 

The techniques of this paper are guided by the principle that they should 
be relevant for non-linear stability problems. For such problems, in the absence 
of symmetry, energy-type estimates have proven the most robust [H]. For La- 
grangian theories like the homogeneous wave equation, such estimates naturally 
arise by contracting suitable vector fields with the energy-momentum tensor 
Taf 3 , to produce a one form Pa- (See |21 for a general discussion.) The diver¬ 
gence theorem relates the spacetime integral of V“Pq, with suitable boundary 
terms. The method can be used to estimate the spacetime integral from the 
boundary terms, but also the future boundary terms from the past boundary 
and the spacetime integral. As we shall see below, both implications will be 
used here. 

Let us recall the situation for the wave equation © on Minkowski space. 
The technique described here was introduced by Morawetz, see ini and Klain- 
erman M By applying the method of the previous paragraph to the Killing 
vector field one obtains the usual energy conservation. By applying the 
method to the conformally Killing “Morawetz” vector field 


K = v^ 


dv 


du 


in the region {1 < t < t^}, one obtains an identity relating a spacetime integral 
in this region to boundary integrals of weighted energy densities on {t = 1} 
and {< = ti}. The spacetime integral can be completely removed by a second 
application of the divergence theorem, which yields additional positive quantities 
on the boundary hypersurfaces. In particular, the weights are sufficient to derive 
decay of the form 0 , assuming that the initial boundary integral is bounded. 
Pointwise decay estimates of the form © can then be obtained by Sobolev 
inequalities, after commuting the equation with angular momentum operators 

n. 

Turning to equation © on the Schwarzschild exterior, we have again a time¬ 
like Killing vector field and thus one immediately obtains conservation of 

®For even the weakest results of |7|, one still requires the analogue of <^ < along 

the event horizon, for an e > 0. 
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the associated energy. Applying as before the above method to the vector field® 
K in the region {1 < t < we again obtain boundary hypersurface integrals 
with a sign, and strong weights, in particular controlling times the energy 
density on {t = ti} D {?; > tt} n {u > tt}. 



On the other hand, even after an additional integration by parts, the spacetime 
term arising from V“Pa no longer vanishes. In the regions r < tq and r > R, 
for certain constants rg, R, however, this spacetime term has a good sign, and 
can be ignored. On the other hand, in the region rg < r < P, it turns out the 
spacetime integral arising from K can be controlled by t times the spacetime 
integrals-summed-arising from vector helds X( of the form 

for a carefully chosen function /f, applied to each spherical harmonic (j)i, (see 
also 12] in connection with vector fields X). In fact, for this bound, the vector 
helds Xi must also be applied to angular derivatives of (j). This leads to loss 
of derivatives in the argument.^ The boundary hypersurface integrals arising 
from Xi are controlled in turn by the total energy. Putting the information 
together from these two vector helds immediately yields energy decay 0), but 
with power —1 in place of —2. 

The boundary integrals arising from the sum of the Xi identities, when 
suitably localized to the future and past boundaries of a dyadic characteristic 
rectangle {tg < t < l-ltg} H {rg < r < P}, can in fact be bounded by tg ^ 
times the boundary integrals arising from K on constant t = to and t = l.ltg 
hypersurfaces. Using this fact, we can iterate the procedure, to obtain that the 
boundary integrals arising from K are in fact bounded, and thus that 0 holds 
as stated. 

The above methods do not give good control near the horizon For 

this we need another estimate, which has no analogue in Minkowski space. This 
estimate arises from applying a vector field of the form 


_ 1 — fx du 

^defined with respect to suitably normalized Eddington-Finkelstein advanced and retarded 
coordinates v and u 

^The necessity of taking angular derivatives is related to the presence of the so-called 
photosphere at r = 3M. 

®This is not surprising, as K is normalized to 
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in the characteristic rectangle depicted. The v length of this rectangle is chosen 
to be of the order of ti. Note that Y extends regularly to the horizon The 
boundary terms arising from Y are related to the energy that would be observed 
by a local observer crossing the event horizon.® The associated spacetime in¬ 
tegral contains a term with a good sign, and terms that can be controlled by 
the sum of the spacetime integrals arising from the X^. As a hrst step, one can 
show using the energy identities for ^ and and the nature of the spacetime 
integral arising from K, that the boundary terms arising from Y are uniformly 
bounded. Then one can go back and, using a pigeonhole argument, extract 
from the term in the spacetime integral with a good sign a constant-u slice of 
the rectangle such that the integral of the energy density as measured by a local 
observer is bounded by Finally, one applies again the energy identity of Y 
in a subrectangle to obtain that the integral of the energy density measured by 
a local observer is bounded by on the segment v = ti depicted. One repeats 
the procedure to obtain decay of 

Note that the procedure above for extracting decay for energy as measured 
by a local observer by means of a pigeonhole argument applied to the spacetime 
integral arising from the energy identity for Y is the analytic manifestation in 
our technique of the redshift effect described previously. Here, the estimate is 
to be compared with the geometric optics argument of the first diagram, but 
where the observers A and B both cross the event horizon, with B at advanced 
time later than A. 

As in Minkowski space, one obtains pointwise estimates (01 from the A® es¬ 
timates via the Sobolev inequality, after commuting the equation with angular 
momentum operators. One should note however, that the necessary L® esti¬ 
mates near the horizon for the angular derivatives arise from the decay rate of 

for the energy flux measured by Y. The boundary terms of K and ^ do 
not contain angular derivatives on U TL~. 

Our use of the vector field Y is inspired by estimates in 0. In that paper, 
one could continue the iteration further to obtain better decay rates than 0 , 
©■ Here on the other hand, the weights on the boundary term arising from 
the Morawetz identity give an upper bound to the amount of decay that can 
be extracted. It would be interesting to explore whether one can surpass this 
barrier using additional techniques. 

A final interesting aspect of our argument is that it does not require in¬ 
verting the Laplacian on initial data or appealing to the discrete isometries of 
Schwarzschild. (In particular, the results here yield an independent proof of the 
classical uniform boundedness theorem m of Kay and Wald.) It may be use¬ 
ful for non-linear applications to avoid techniques so heavily dependent on the 
exact staticity. In this sense, the argument given here is perhaps more robust. 

®More precisely, exactly the part of the energy not seen by 
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2 Schwarzschild 


Let (M^g) denote the maximally extended Schwarzschild spacetime with pa¬ 
rameter M > 0. This manifold is spherically symmetric, i.e. the group SO{3) 
acts by isometry. We recally briefly the usual description of the global causal 
structure of M, via its so-called Penrose diagram. (We refer the reader to 
standard references, for instance nao 

For spherically symmetric spacetimes, recall that Penrose diagrams are just 
the image of global bounded null coordinate systems on the Lorentzian quotient 
Q = A4/SO{3) viewed as maps in the obvious way from Q In the case 

of Schwarzschild, the Penrose diagram is depicted below: 


B 



The curve S depicts the projection to Q of a particular choice of complete 
Cauchy surface S C . We will call the sets J~ ) C ) and J~ {1 ^) C 

the exterior regions.^® We call the future event horizon correspond¬ 
ing to the end and Hg the future event horizon corresponding to . We 
could also consider Tig as the past event horizon corresponding to , and thus 
may denote it alternatively by . 

Defining a function r : Q ^ R 

r{q) = y^Area{q)/4:Tr, (6) 


the two exteriors are each covered by a coordinate system (r, t) so that the 
metric g may be written: 



dt^+ll 



-1 

-I- r^dCT§2, 


(7) 


where dagi denotes the standard metric on the unit sphere. The vector field 
^ is clearly timelike Killing. These coordinates break down on U H^. The 
Schwarzschild solution was originally understood as the spacetime described 
by the expression Q. As explained in the introduction, the realization that 
this was actually just one of the exterior regions of a larger spacetime, took a 
surprisingly long timed^ 

^^For an explanation of the notation etc., and more about Penrose diagrams, see 

the appendix of 1^ . 

^^Note that the expression Q also describes the metric in the black hole interior fl 

\ (^B ^ ^ spacelike Killing. 
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The sphere is known as the bifurcation sphere of the event horizon. 

If we chose S to contain n T-L\, then 

S' = 

is a Cauchy surface for the exterior region J~ il^) n J+(X^). In particular, 
solutions of the linear wave equations on J~ {X\) n ) are determined by 

their data on S'. 

The above means that if one is only interested in the behaviour of solutions 
to lO in J~ il\) n J+ (X^), one can study the problem in any coordinate system 
defined globally in this region, in particular, so-called Schwarzschild coordinates 
(r,t). For convenience, we shall in fact use a null coordinate system {u,v), 
which “sends” T-l\ to u = oo and Hg to v = —oo. These coordinates, so-called 
Eddington-Finkelstein retarded and advanced coordinates, will be described in 
the next section. 

The reader should not think, however, that imposing such a coordinate sys¬ 
tem removes the geometry from this problem. For one must not forget that 
the correct assumptions on (j) are those expressable geometrically, not those 
that happen to look natural in the chosen coordinate system. When written 
in local coordinates, these assumptions would be difficult to motivate without 
knowing the origin of the problem. Moreover, the same can be said for the 
form of many of the techniques used here. In particular, the form of the vector 
field Y of Section 0 is best understood by passing to a new regular coordinate 
system on Ti^. In the coordinate system chosen here, Y appears asymptotically 
singular. 

With these warnings in place, we turn to a description of two related coor¬ 
dinate systems covering J~ {X\) n J+(X^). 

3 Eddington-Finkelstein and Regge-Wheeler co¬ 
ordinates 

Let (r, t) denote the Schwarzschild coordinates of 0. Define first the so-called 
Regge-Wheeler tortoise coordinate r* by^^ 

r* =r + 2M log(r- - 2M) - 3M - 2M log M, 

and define retarded and advanced Eddington-Finkelstein coordinates u and u, 
respectively, by 

t = V + u 

and 

r* = V — u. 

^^For instance, one should not restrict to (p of compact support on S', for this would mean 
that (j) necessarily vanishes at n 'H\, an actual sphere in the spacetime. See also the 
remarks in the Introduction. 

Coordinates are together known as Regge-Wheeler coordinates. We have centred 

r* so as for r* = 0 to correspond to r = 3M. This is the so-called photosphere. 
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These coordinates turn out to be null: Setting fj, = 2M ^ the metric has the form 

—4(1 — fi)dudv + r^da^. 

We shall move freely between the two coordinate systems {r*,t) and {u,v) 
in this paper. Note that in either, J~{X\) n is covered by (— 00 , 00 ) x 

(— 00 , 00 ). By appropriately rescaling u and v to have finite range, one can 
construct coordinates which are in fact regular on and By a slight 

abuse of language, one can parametrise the future and past event horizons in 
our present {u,v) coordinate systems as = {( 00 , ?;)}„g[_oo,oo) and = 

{(^; ^^)}lie( —00,00] ■ 

Finally, we collect various formulas for future reference: 


2M 



9uv = (ff™)-' = -2(1 - m), 
dyT = (1 - m ), duT = -(1 - n) 
dt = dv + du, dr* = dv — du, 

d 1 / d 9 \ 
dr* 2 \ clu du) ' 

dVolM = r^(l — fi) dudvd(Tg2, 

dVolt=const = r'^ y/1 - fj, dr* dcrs 2 , 

= -(1 - {d^^ - r-2a,.(r2a,.V')) + y^yA<P- 

Here y denotes the induced covariant derivative on the group orbit spheres. 


4 The class of solutions 

Let S' denote the surface {t = 1} say in the exterior, and let 5 be a Cauchy 
surface for M. such that S C J“(X+) n J+(X“) = S'. Let N denote the future- 
directed unit normal to S. 

We proceed to describe the solutions cj) : A4 ^ M. oi the wave equation 

= 0 ( 8 ) 

on A4, which we shall consider in this paper. Given (j), N(j) on S, define the 
quantities 

Eo = ^ f /r2(l-^)-5 (^(at(r*y»)2-f (9^.(r*yV))2 

+(1 - {^,r*,<7s2)dr*das2, (9) 
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2 

£’i = ^ / f r'^ (u^{duir^f"(l3)f +v^{dy{r^Y(l)))‘^ 

i=0-{o £2 

+(1 - fj,){u^ + ?;^)|y(r*yV)P) (1, r*,u;)dr*das ,2 


E 

i=0. 


rM(at(r*yV))^ + (5,.(r*yV)f 


£^2 = 


+ (1 -/^)iy(»'*yv)p) {l,r*,u;)dr*das 2 

4 oo 


+(l-M)(u"+7;")|y(rY 

5 




(1, r*,u;)(ir*(icrg2 


E 

2 = 0 


(54r*y»)2 + (a,.(r*y»)2 


^-oo S2 

+(1 - d-)\f{ry''4>)\‘^] il,r*,uj)dr*dcrs^ 
2 


E 

2 = 0 . 


r2(l - {dt{ry^c^)f + (9,. (r*y»)^ 


+(1 - (l,'r*,a;) • dr*das2. 


( 10 ) 


( 11 ) 


Our weakest result, namely, the uniform boundedness, requires the boundedness 
of Eq. Our energy decay result will require the boundedness of Ei and our full 
pointwise decay results will require the boundedness of i? 2 - 

We note that the boundedness of Eq follows from the statement that cj) is C^, 
Ncj) is on S, and that (p and decays suitably at spacelike infinity, for instance, 
if (j) vanishes identically in a neighborhood of Similarly, the boundedness of 
El follows for suitably decaying p and Np, and finally the boundedness 
of E 2 follows from C® </>, etc. There is no assumption of the vanishing of </> on 

H+nn-. 

We will state a coordinate version of Theorem o in Section 0 As this 
theorem will refer to the energy flux defined by the Killing vectorfield we 
will first need some general results regarding conservation laws. These will be 
given in the next two sections. 


5 Conservation laws 

As discussed in the introduction, the results of this paper will rely on estimates of 
energy-type. Such estimates arise naturally in view of the Lagrangian structure 
of the theory. We review briefly here. 
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In general coordinates, the energy-momentum tensor for (j) is given by 
Tap = da<j)dp(j)- ^gapg'^^d^<l)ds(t). 

This is divergence free, i.e. we have 

= 0. (12) 

For the null coordinates we have dehned, we compute the components 

Tuu = 

Tvv = 

Tuv = = (1 - M)iy<(>r 

Let denote an arbitrary vector held. Let n“ denote the normal to a hyper- 
surface^^. Let TTy^ denote the deformation tensor of V, i.e., 

7r“^ = i(V“V'^-hV^y“). (13) 

We will denote in what follows the tensor iVy^ just by 7r“^. In local coordinates 

we have the following expression: 

7h/37r“^ = {{du(l))'^dv{Vv{l - g)~^) + {dv(j)fdu{Vu{l- 

+m\'^{duV, + a„K)) - :^(K - K)(iy0p - rM- 

2r 

Let 5 be a region bounded to the future and past by two hypersurfaces Si 
and So, respectively. Let P“ = g^'^TpsX^. The divergence theorem together 
with m and m gives 

f TapTT°‘'^dVols = / gapP^n^dVol^^- f gapP'^n'^dVol^g. (14) 

t/ 2i ^ Sq 

6 Conservation of energy 

Let us apply the above to the vector held As this is Killing, by and 
m, the associated vector held P is divergence free, so there is no space-time 

^"^Note the convention in Lorentzian geometry 
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term. The (negative of the) boundary terms on constant t-hypersurfaces are: 

OO 

= J J+ ^Tuu + ^Tuv'^ ■ r'^^/T^^dr*das2 


— OO 5 
OO 


VI - M 
1 

4VW 


{{dycl))^ + + 2(1 - n)\f(l3\‘^) ■ r^\J\ - fi dr* da §2 


— OO « 
OO 


/ / + + {1-■r'^V^^^dr*das2 


on constant-M hypersurfaces are given by: 

^1)2 p 2 

I VI Js'^ 4(1 — /r) 

_ n f 1 

Jvi J§^ 4(1 - fj.) 
and on constant-u hypersurfaces are given by: 

1 


pV2 P 

F^{{u} X [ui,r; 2 ]) = / / 777 -— s{Tvv +Tuy) ■ r^(l - fi)dvdas 2 

{{dy(j)f + (1 - /r)|y(/)n • r^(l - 7t)dwdcrs2 


Fcj>{[ui,U 2 ] X {v}) = f f —-(r„„ + T™) • r^(l - 7 t)dM(icrs 2 

Ju, Js^ 4(1 - /rj 


Jui Js^ 4(1 - n) 

The identity Cl applied to the above shows that 

F^{to) = F^{1) = Etj) 

and that the fluxes satisfy 

F^{[ui,U 2 \ X {v}) < F^, 
F^{{u} X [ui,ii2]) < F^. 


{{du(l)f + (1 - d)\f(l)?) ■ ?'^(1 - /r)dudcrg2. 


In particular, we can define a function on Q n J“ (X+) n J+ (X~ ) by 

m^{u,v) = F^{[u,oo] X {u}) + iV({oo} x b^oo])- (15) 

It is clear that the bound 




(16) 


follows immediately. 
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7 Coordinate version of the main theorem 


Theorem 7.1. Let ^uj), be functions such that the quantity Eq 

o/ 0 is bounded}^, and let (p be the unique solution o/0 on J~ (T^) O {T~) 
with (p(r*,l,uj) = (j)o{r*,uj), dtpir*,l,uj) = (pi(r*,uj). Let vu^ be as defined in 
Then there exists a universal constant C such that 


\(j)\ < CEqU 2 . 

If the quantity Ei of GB is hounded, then 

070 < CEi{vf'^ + 

If the quantity E 2 of m is bounded, then 

\cP\ < CE2vf} 


for all {u,v,iij), while 

\rfi\ < Cf,E 2 (l + |u|)-^ 


(17) 


(18) 


(19) 

( 20 ) 


for all r > R > 2M, t>l. 

In view of our previous remarks, it is clear that Theorem o implies Theo¬ 
rem 


Remark 7.1. The number of derivatives required in the definitions of Ei,E 2 
can be reduced by a slight refinement of the analysis in Section\^ We shall not 
pursue this here. 


8 The vector fields Xi 


In this section we shall define, for each spherical harmonic (p^, a vector field Xi 

by 




--fg— + -h— 

r^du 



( 21 ) 


for some function fi = fi{r*) to be determined later. We shall show that the 
function fi can be chosen so as for the spacetime term corresponding to Xt 
to be positive, and so as for the boundary terms arising to be controlled by 
the usual (conserved) ^-energy, after application of a Hardy inequality. (See 
Proposition irr^ from the next section.) These V^-energy identites can then 
be summed so as to yield an identity relating a positive spacetime integral and 
boundary terms controlled by the usual ^-energy. (See Section 

the obvious sense, i.e. when (p is replaced by (po and dt4> is replaced by (p\ when 
evaluating m 
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8.1 Identities 


We first collect various identities for vector fields of the form m for a function 
/. For functions of r*, let ' here denote The spacetime integral given by 
the left hand side of da in a region TZ = {to < t < ti} is: 





( 22 ) 


while the boundary terms are given on a constant t-hypersurfaces by: 

OO 

J J - fidr^, das2. 

-OO §2 


We have the identity 

if{n)=Ef{t,)-Ef{to). (23) 


To produce an identity with terms for which we can control the signs, we wish 
to replace the spacetime integral with a new integral obtained by applying 
Green’s theorem to its last term. Defining 








( + 
V2(1-m) 






{1 — fj,) dt dr* da§ 2 , 


/ 


(24) 


we have by Green’s theorem that 

CXD 

J J + 2^—^'^ dt(l)(l){ti,r*,uj)r'^ dr* das 2 

— OO 
OO 

- J j ^ (^f + 2 ^ ^ dt(l)(t){to,r*,uj)r‘^ dr* das 2 , 

-OO §2 


and thus we have the identity 

I^{TZ)=E^{h)-Ef{to) 

where 

OO 

Ef{U) = E^{U)+ J J + 2^—^'^ ^t(|)(|){U,r*,UJ)r'^^/^E^J{dr* das2. 

— OO 

(25) 
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Finally, let us compute the expression with the d’Alambertian on the right 
hand side of (EH). Since, for a function ij) dependent only on r, we have Otj; = 
(1 — it follows that 


°(/' + 2 V^/) = 


1 2 2 
+ -/" + -r 

1 — /r r r 


1 - /i 


1 - M 


r 


+ 4(1 - m)/' - 4(1 - - 4^'/- 


In view of the relations 


(1-/^)^ 


and 


I-A” I o (l-/^V 

,^2 


we obtain the expression 


□(/' + 2^/) = 


1 p/ll , ^ J-ll _ 'IM pI 

1 — fi r(l — y) 


(l-/r)r 




- A /• (26) 


8.2 The choice of fi for 1>1 

In this section, we shall choose fi for the higher spherical harmonics i>\. 

Our goal is to make the spacetime integral (TZ) positive, and the bound¬ 
ary terms controllable by the energy Thus, one would like the function 
to be bounded, and, in view of the first term in lEH), should be positive. 
For this to be the case, however, /"' must become positive in a neighborhood of 
the horizon. Examining, however, EH) and EH - one sees that this term enters 
in with the wrong sign, and, what is more, it is multiplied by (1 — 
and thus seems to dominate the only term with the correct sign to cancel it, 
namely f^. 

To generate a term which can indeed cancel this “bad term”, we must borrow 
from the term in EH with dr* 4>- This is possible as follows: For any function 
/3 = (3{r*) decaying to 0 as r* ^ ±oo, we note the following identity: 



— {dr*Ar‘^{l - A)dr* = 
I- 



/3 + 


— {dr*(l) + - fi)dr* 

l-/i 

-r^—d' + -f'p) - Adr*. 

l-/i r J 
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Let us set M = 1 (i.e. = ^) and, for a sufficiently large constant a to be 

determined below, let x denote the coordinate 

X = r* — a — 1. 


Define 


/3 = 


so that 


1 — /J, X 






^ r 1 + x"^ {a^ + x^)"^ 


/3' — /3^ + -(1 — fj,)P — ——— — 

r r (a^ + x‘^Y 


We may thus write 



(4/r - 3) r"^ {1 - fi) dt dr* das 2 . 


For each spherical harmonic number in view of the relation 
[ |y(/)^pr^dCTs2 = f + ^\ cj)e)'^r'^das2, 

we may rewrite 



+ (1 “ d') dtdr* das 2 . (27) 

The expression multiplying fi(p'j in the last term above vanishes at a unique 
value of r*. Denote this by —j£. Note that 


-4/3 < r*{8/3) < < r*{3) = 0. 

We may now define fi = fi{r*) by setting 

fehim = 0 , 
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{feY = {c? + x^) . 

Let us drop the I in what follows. We have 

‘lx 


r = - 


+ x'^y 


Denoting by 


we compute 


F = - 


r = 


1 




\2 ■ 


(a^+a;^) (a^ + a;^)" 

4/"a; 4/'a2 


4(1-m) 


F = 


r + 


o? + x^ (o? + x^y 


2 2 
— a 


2(1 m) (^a^ _|_ j.2'^ 


3 ■ 


Note that F is nonnegative for x < —a, and x > a. 

We would like to show that for a sufficiently large a, independent of we 
can dominate the second term in pointwise by the last term, i.e. we want 
to show 

-F< 


2r3 


/• 


(28) 


In view of the sign of F, it suffices to consider —a<x<a. 
For —a<x<a then, we estimate 


/ = 


r f > 

J _Q;_ l_-y 


X + a Fl 

2a2 + 2Q; + r 


To show PH|i. it suffices then to show 


1 ) 


1 — 2)^ 3 — 4/r\ a: + a + l 


> 


1 


2 2 
a — 


(29) 


2r3 J 2a^ + 2a + 1 2(1 —/r) (a^+a;^)^ 

Consider first the region —a < x < —laji. For a sufficiently large, we have 


1 


a^ — a;^ ^ i^/l a + x 


2(1-m) (a2 + a;2)^ ' 2 ' 

On the other hand, since for x < —laj?) we have, say r < 5a/12, and since 




we have (2 — 3/i)(r*) > 1/5 for r* > 1. Then for £ > 1, 

«(f+i)T4£^_4±£±i_>* + “ 
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for sufficiently large a. This gives (P|l for the region considered. 

Consider now the region —2a/3 < x < a. Note that as a ^ oo, we have 
r~a; + Q;, /i~0in this region. Thus, we have 


£{i+l) 


2 —3^a; + a+l + 7 £[£+l) x + i 


while 


2j- 3 2a'^ + 2a + 1 2 (cc + a)^a^ ’ 


1 (a — x)(a + x) 


2(1-m) {a^ + x^y 2 (a2+2.2^3 

To show (HU for sufficiently large choice of a, it suffices then to show the bound 

(a — x){x + a)3a^ 


2(a^ + x"^)^ 

For X < 0, 1301) is immediate from 


< 1 . 


{a-x)a^ ^ y/2{a^+x'^)a^ ^ 


2(a^+x^)3 2(a^+x^)3 

For X > 0, we have on the one hand 

(a — x)(x + a)3a^ ^ (x + a)^a^ ^ 1 f ax + 


2(a^ + x^)3 2(a^ + x^)3 2 \ x^ + a^ 

On the other hand, 

1 


ax + < ^x^ + ( 1 + 


2q 


a\ 


Set 


Then q = 1 + and 


q = 2 + -, q-2q-l = 0. 
9 


■ + a^ < i(l + V2)(a^ + x^). 


The bound 1201 then follows from the inequality 

((1/2)(1 + V2))'<2. 


(30) 


8.3 The case ^ = 0 

For .^ = 0, we have the identity 

2Th^7r“^ = {{dtfjiof + [dr*4)of) + ^ {{dt4>Qf - [dr*(j)Qf‘) . 
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Applying 12311 with the sharp cut-off function /q = X{-oo,r*)j so that /q = 
—6{r* — R*), we obtain 





to §2 


((5t<^o)^ + {.dr*4>of) {t,R*,L0)r'^{l - fi)das2dt 


-b 


II 


to —oo S2 


{dt4>of) »'^(1 - fi)das2dr*dt < 


Let us denote 


ti r* 


F{r*) = J J J- n)das 2 dpdt. 


Then, in particular, 


to —oo 


r 1 


Hence, 


p(* 


2Y^F\r*)<Fir*)+4E^,. 


Fir*)) <4e 




2r (1 p,)Erp^. 


Note that fl^2r ^(1 — p)dr* = log/r Also observe that F(—oo) = 0. 
Therefore, integrating we obtain 


and thus, 

It now follows that 

tl 


p'^Fir*) < 4E^, [ fi^df < AE^^, 

J2M 

Fir*)<Ap-^E^,. 


J J 1 ft,r*,a;)r^(l - p)dadt < Ap ^E^^. 


to S'2 

Multiplying by and averaging over r* we obtain 


f^O J_ 
''do 




to — oo §2 


(1_^)(1+|^*|)1 + 


({dt4>)^ + (9r-*0)^) ?"^(1 — p)dadr*dt 


(31) 
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8.4 The identity, summed 

We introduce the notation 


E 


X>i 

0 




and 

9 / > Vt 

i>l 

Our final energy identity for the totality of vector fields Xi, I > 1, can be 
summarized by the statement 

0 < - E^>-^{to). (32) 

Defining 

=^"+d“‘. 

< = (2+)E*. + <-‘, 

we have the inequality 

lfiTZ)<\E^{h)\ + \Efito)\. 


9 The vector field Y 

In this section, we shall introduce a vector field which will give good control on 
the solution near the event horizon Although the computations are done in 
Eddington-Finkelstein coordinates, the reader is encouraged to compare them 
with computations in a null coordinate system which is regular on the event 
horizon. 

We set Y to be the vector field 


Y = --f— - . 

2^ du 2"^ dv 


where 


/ = 


a(r*) 

1 - Ai’ 


f = P{r*), 


for functions a and f3 do be defined later. Below ' will denote We have 


rj. ^5 fafi A {dy(/))^ , 


4(1 

-\\n? 


4(1-m)' 

a' (/3(1-Ai))'\ 1 


1 - M 


1 - M 


a 


2r \1 — 


- P dupdyp. 
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Integrating in the characteristic rectangle TV = [ui,U 2 ] x [v\,V 2 ]^ we obtain 
^^({“ 2 } X [vi,V2]) + {[ui,oo] X {^ 2 }) 

= {T^') + fJ {ui X [vuV2]) + {[ui,oo] X {wi}) (33) 


where 


V2 

F^i{u} X [vi,V 2 ]) = j J + P{dy(l))'^){oo,v)r'^dvd(T§ 2 , 
S2 


Fjiiv} X [ui,U2]) = J J + {l-r^duda§2, 

§2 


'^2 U2 2 

{F') = 111 (^- ( 1 ^ -«') + +\n\^ («' -m- m))')) 

Vl S2 

2 \ 2 

+ -(a — /3(1 — /x)) 9„(() dy(l) r du dv da§i. 
r ) 

We define a, (3 as follows. Let tq > 2M be a constant sufficiently close to 
2M, to be determined by various restrictions that follow. Set a = 1 and /3 = 0 
on the event horizon. Furthermore, require that a, f3 both be non-negative 
functions supported in the region r < 1.2ro, with 

a'ir*)^C{l + \r*\)-^-, f3{r*) ^ C{1 + \r*\)-^-, Vr < ro 

for some constant C > 0. 

We will always apply the above identity in [mi,oo] x [r'i,n 2 ]. We have then 
that 

X [x;i,'(;2]) = J J\f(l)\‘^{oo,v)r‘^dvdas2. 

VI §2 

Let us set 

V2 00 

iJ{n')=J J J ~ ^') + {dvpfp' + {a' - (/3(1 - r‘^dudvdas 2 . 

VI 111 §2 

The quantity tq will be chosen sufficiently small so that all terms in the above 
integrand for ij are nonnegative in the region r < tq, and so that moreover 

( 34 ) 
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8r-^(l-^) </3' (35) 

in r < tq. In addition, we will require of tq that 1.2ro < 3M (see Proposi¬ 
tion and (151) . 

Defining 

V2 OO 

lJ(TZ')= j J J ^{a-P{1-^i))du4>dy4)r'^ dudvdas2, 

VI Ul §2 

we may rewrite (Pl) as 

Fj ({oo} X h, U2]) + Fj ([z.1,00] X {U2}) + il [F!) 

= X [vi,V 2 ]) + FJ{[ui,oo] X {ui}). (36) 

In our argument (see Proposition II 1.51) . the terms and in the region 
r > tq will be controlled by . The relation of F^ with the integrand of 
(as exploited in Proposition 1 11. (ill can be thought of as the manifestation of the 
red-shift effect as measured by two local observers. It is to be compared with 
the fact that, in equation (51) of |5] for the factor in front of ^ on the right 
hand side is bounded above and below away from 0 uniformly in the region 
r < ro- 


10 The vector field K 


In this section, we shall define the Morawetz vector field K. 
We set K to be the vector field 


We have 


^ _ _1 2 d _ 1 2^ 

2^ du 2^ dv' 


K- = if" = if„ = (1 - tJ.)v^Ky = (1 - 


We compute 




4r 2r 1 ^ ' 


47- 


Let TZ = {to < t < ti}. The identity itTl gives 


I^iTZ)=Efito)-Efih) 


where 


/f(7f) = [ T^pn^^dVols 

Jn 




'to 'J —oo ' 


2 4r 


2r 


Ar } 
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and 


= / gajsn^'P^dVol^t^t^} 

1 


2vr^ 

• ^ dr* dtTg 2 


{K^T,u + K'^Tuu + K^Tu. + K'’T,,){r* ,U,uj) 


1 


VI - M 


{U^TuU + V^Tyy + (it^ + V^)Tuv) 


yjl — ^ dr* das 2 
1 


/—oo</S^ d-x/T J1 
Y^l — /z dr* 


(u^idu^y+v^d.^y + (1 - m) V + vymn 


In view of the identity 

f 'il)U[(y)dVol-ji — [ {Otp)cj)'^dVolTi 
Jtz Jn 

= [ {2g°'f^cj)\7a(l)n0'ilj - g°‘^Waipnf3(y)dVol{t=ti} 


l{t=to} 


(2g“^(^Va V/3'*/' - g°‘^'^ayni3(y)dVoyt=to}, 


applied to 

we obtain the identity 


y = 


tr*{l-g) 
4r ’ 


iy{ji) = Ey{ty-Ey{ty 


(37) 


where 


pti pOO 




'to j —oo 


+VVIV 2 + 


tm\ 

r*(4/r- 3) 


2 V _ 

2 4r 


and 

^Kf 


Ef{U)=E^{U)+ -y= 


1 /tr*(l-/i) 


M 


2r 


2r 


■ r^ (1 — pi) dt dr* dfT §2 (38) 


(/)dty - Vl - M c?crs2 ■ 


For the analogue of the following argument in Minkowski space, see El Let 
S = vdv + udu and 5 = vdv — udu- Then 

v'^{dy(j)y + u^iduy)"^ = i {{sy)'^ + {sy)"^), 
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S = 2{tdt+r*dr*), S = 2{tdr^+r*dt), 

t 

2tdt(l) (j) = S(j) — 2r*dr»4>, 2tdt4> = ~S. — 2— 

Thus, by integration by parts, we obtain^®: 

J 2t ^^~ dt<i)(l)r‘^dr* = J ~ + 8^* ((1 -/r)r(r*)^) dr* 


(1 - /i)r^ ( —S84>+ ( + 2— ) (j)^ ) dr*. 


It now follows that 


f + 2t— — dr*d(j §2 

Js 2 V4 r r J 


(1 - {{IS4+ dr*das2. 


We may thus write 


E^u) = I ^^^(\{i+fi){s4? + As4r+a- 


VI - M V4 




+ (1 - /r) (^2'^'^ ~'^j + j ■ d(TS2.{m) 

In particular, E^{ti) is nonnegative. 

11 Comparison estimates 

Let us introduce one hnal spacetime integral quantity: For X a region denote 

i^[x)= [ {\dr^4\'^ + \dty + {i-fim\'^)dVoi. 

Jx 

We have the following 

Proposition 11.1. If X is a rectangle defined by 

X = {to <t <ti}r] {2M <ro <r < R}, 

then 


I^{X)<C{iyx)+ll{X)). 
Moreover, for either R < 3M or vq > 3M, 

I^{X) < ClfiX). 


(40) 


(41) 


^®Note that there is no problem near the bifurcate sphere in the integration by parts argu¬ 
ment due to the presence of the exponentially decaying (1 — fi) factor. 
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Proof. From the definition of the functions /^, we have that for ^ > 1, 
r ^2 

X, 


Jx (1-m) + 

and thus, summing over £ > 1, we obtain 


rdVol < Cll{X) 


^>1 


,dVol < 


Jx (1 - /r) (1 + {r*r) 

Now using this bound, we can apply again with the function 

dp 






(42) 


(43) 


(44) 


Note that this choice ensures the non-negativity of the space-time integral con¬ 
taining the angular derivatives. We obtain 


\dr>cj) 


>1 


(l-^)(l+(r*)2) 


dV ol 


f{r*){2-3pm 


>1 


lx 


2r 


-dVol < 


Applying (1^ . using the representation 


(45) 


° + (1 - M)iy<(>>in (46) 

and the function f{r*) = f{r*) ■ (2 — 3p){2r)~^, we obtain from (1^ and 
the bound 

[ /(r*)(2-3/i)(2r)-i|5t<(.>ipdroZ<C/^^^(A). (47) 

Ja: 


On the other hand, by m, we clearly have 




lx {1 - p) (1 + {r*r)^+ 


dV ol 




lx (l-/r)(l + (r*)2)i+ 


dVol < Cll{X). 


Estimates (1^ . (in and m immediately yield 


(48) 


f /(r*)(2 - 3p)\dt(l)\'^ p'^ 
lx 2r{l + {r*ry+ 


dV ol 


\dr*(t>\'^ p'^ 


lx (1 - ^) (1 -I- (r*)2)^+ 


dV ol 


+ 


f nni2-3p)\n\\ yoi<cix 
Jx 2r - ^ 


(49) 


This clearly gives 6D- 

Applying our argument again to the solution = flp of the wave equation 
= 0, where f2 is an angular momentum operator, in view of the fact that 
this solution has vanishing zeroth spherical harmonic, we obtain as in m 




lx (1 - p){i + {r*ry 


rdVol < ClfJX). 


(50) 
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Here actually denotes the sum of I applied to fl(j) where H range over an 
appropriate basis of angular momentum operators. Then, using and (EHl), 
and arguing as before we also obtain that 




(l-/r)(l + (r*)2) 


rdVol <C{lfiX) + 


The estimates (P)l. PHl and iu together give roni . 
Proposition 11.2. Consider the hypersurface {t = ti}. 


'—I 

(51) 

We have 

□ 


\E^iU)\<CE^. 


Proof. In view of the Cauchy-Schwarz inequality applied to the boundary terms 
defined by (123, it suffices to to establish the Hardy inequality 

OO OO 

j y(l + |r*|)“^|(^|^r^(ir*dcrs 2 < C J j\dr» r^ dr* da§i. (52) 

— OO §2 _oo g2 

Define the function f{r*) by solving the equation 


^'“iTTF 

with the boundary condition /(—oo) = 0, i.e., 




Here ' denotes Clearly. 

n.r * 'J • 


fin 


\ (for r* —oo) and /(?'*) 


(for 


+oo). 
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Now, 




= -2 


< 2 


-OO 

CO 

-OO 

ifir*)f 


= 2 


dr* 

f'{r*) 

-OO 

OO 

I J^J’{r*M\^dr* 

— CO / 

{{l + \r*\)f{r*)f 


§2 


\dr* (l)\'^ dr* 


OO ' 

I y^(l + |r*|)"2|(/)pr2dr* 


Inequality follows immediately from 

{a + \r*\)f{r*)f 


< Cr\ 


□ 


We will denote by the total energy associated with the solutions of 
the wave equation 4>^ = fl(p, obtained by applying a basis angular momentum 
operators ft to and summing. We will also set = E^ + E^^ etc. 

Proposition 11.3. Consider the rectangle X of Provosition \mi We have 


I 4 X) < 

If X C X {r = 3M}) = 0 then 

U(X) < CE^. 

Proof. This follows immediately from the previous two Propositions. □ 
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Proposition 11.4. Consider the rectangle TZ = [mi, oo] x [vi, V 2 ], and the shaded 
triangle TZ depicted. ° 



We have 


Pj {{Ul} X [VI,V2]) < C{ZU^{UI,V2) - 

w^{oo,V2) - W^{u{v2),V2) < CF^ ([uijOO) X {?;2}), 
00 



das,2du < CFJ ([mi,(X)) x {i;2}) 


U{V2) 


v{u) 

J J \f(l>\'^das 2 dv < F^{{u} x [i;!,«(«)]), 

1^1 S2 


Vu > ui 


where u{v 2 ) and v{u) are defined so that r{u{v 2 ),V 2 ) = rg and r{u,v(u)) = ro 
respectively. 

Proof. This follows immediately from the definition of Y in Section O □ 


The next proposition shows that in the identity (El generated by vector 
field Y, the space-time terms without a sign can be controlled with the help of 


Proposition 11.5. With TZ, TZ as above, we have 

F^(.{oo} X {vi,V2}) + FJ{[ui,oo) X {U 2 }) -f 

< C {l^iTZ) -I- -F^([ui,oo] X {ui}) +w^{ui,V2) - w^{ui,vi)) . 

Proof. By Proposition nn we have that 

X [U1,U2]) < C{w^{ui,V2) - W,j,{ui,Vi)). 

From EH), it would suffice then to show that 

i^{1z)<cifijz) + \il{iz). (53) 

Let US decompose: 

il{TZ)=il{TZ\TZ)+il{TZ). 
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By an application of Cauchy-Schwarz and CT . (1^ . we have in the region TZyR 
that 


^(a - /3(1 - < I 

^ 1 f {ducj))^ fafi 

- j' 

so integrating, we obtain 

i^{n\n)<^i^{n\n). 

On the other hand, by Pronosition ll 1 .H and the condition 1.2ro < 3M, we easily 
see that 

iJ{TZ) < C'(/^(7^)). 

Since 

0<iJ(TZ\TZ)=iJiTZ')-iJiTZ), 

and 

\iJ{TZ)\<Cilf(TZ)), 

follows immediately. □ 

The next proposition can be thought of as a pigeonhole argument, which 
will allow us to pick values of v where the boundary terms generated by the 
vector field Y gain additional decay. 

Proposition 11.6. With TZ as above, we have 

inf FT ([mi, oo] X {u}) < C{v 2 - vi)~^iY{H) + C(w0(mi, u) - w^{u{v), v)). 

Vi<V<V2 ^ ^ 


Proof. This follows immediately from the fact that f > {h — a) inf /, and the 
inequality 


F^([ui,oo] X {u}) < 


r f (iMt 
Jui dS2 \(1 - M) 



+ {dy(j})‘^f3' + (o' - (/3(1 - m))0) f^duda 


+ w^{ui,v) - w^{u{v), v). 


□ 

The next proposition estimates the positive part of (i.e. the terms with 
the wrong sign) by I{X), with a certain loss: 
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Proposition 11.7. Let TZ' and X denote regions TV = {to ^ t < ti}, X = 
TZ' n {to <r< R} depicted below: 


t = ti 



t = to 

where R is sujficiently large. We have 

I^{1Z')<hI^{X). 

Proof. Recall that 


ti OO 


t , 


inn') = 


4 7-2 

We have required of tq that 

Tg < max I — 


3/i 2 

1 H-^- r 


to —OO 

^t (2p ^ (4/r-3)^^^ 


2rn 


2r 

r* \ (fn r"^ {1 — fi) dt dr* (i(Tg 2 . 


2rn 


4^*0 ~ 3 3/ro — 2 


< 0 . 


(54) 


We require of R that R is sufficiently large, to be determined below, and, in 
particular, that R > r 2 where r 2 is the infimum of all r satisfing, for all f > r, 
the inequality 

2f 2f 


r > min — 


4^ — 3 3^ — 2 


> 0 . 


(55) 


Note that r 2 < infty. The choices (EU, (ESJ ensure that in the regions {r < rg}, 
{r > R}, the integrand of is nonpositive. Thus, we have 

if{iz') = i^{iz'\x)xif{x) 

< ifi.x). 

On the other hand, one sees easily that 

(^1 + < Chl^iX). 

It will then suffice to prove the inequality 

f i ( 4 m - 3 )m 

Jx 4: \r^ r3 


r* dVol < -elf in' \ X) XCtil^iX). 
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for some sufficiently small e > 0. The function 

is positive only in a subset of the region of 2M < tq < ri < r < r 2 < R. 
Moreover, H{r*) behaves like r~^ as r —> +oo. The desired result will then 
follow from the following inequality: 


'2 '3 '3 

J j cj)^r^{l—fi)d(7s'^dr*<C j j[(j)'Yr^{\—^i)da§2dr*+c J j r^{\—^)da§2dr* 


S2 


for some small constant c > 0 and some large < R. This in turn easily follows 
from rescaling the following one dimensional estimate: 


Lemma 11.1. 

such that 


For any c > 0 there exist sufficiently large constants C and A 

/ 2 nA pA 

(f^{x)dx<C J {(j)')^dx + 2 c J x~^(j)^dx. 


Proof. We write 


(/>2(x) 


4'^iy) ~ ^ [ 4’'4'dz < 2 

J X 




Thus 

(j)^{x)dx < c~^ J z(j)'^dz + c^ j z~^(j)‘^ dz + (^^{y). 

Dividing by y and integrating in y in the region [2, A\ with A = we obtain 

/ 2 pA pA pA 

(j)^{x)dx < c~^ J zcf'^dz + c J z~^(j>'^ dz + j y~^(j?{y)dy. 

The desired conclusion follows immediately with C ~ Ac~^. □ 

□ 

Proposition 11.8. For the region TZ' as above, we have 
E^{h)<C{I^{n') + E^{to)). 

Proof. This follows immediately from (EH). □ 

Proposition 11.9. On a constant ti hypersurface, we have the bounds 

VI - M 

+ (1 - y){u^ + V)|V(/)p) • - y dr*da ^2 < E^{ti), (56) 
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(57) 



- ndr*das2 < 

\/l — fJ,dr*d(T§2 < CE^{ti). 


(58) 


Proof. The bounds (1^ and follow immediately from 
For , write 

J 2t- -—— dtfxfr'^dr* = J + t'^dr* {{I — ^.t)r) (f^ 

=y(i - m)’’^ {^S4)(j>+ 

=y (1 -m)’’^ 


and use dSl. □ 

We now compare the boundary terms E^(ti) generated by K with the flux 
of ^-energy. 


Proposition 11.10. WithTZ' as before, let {rl,ti), (r*,ti), be so asforti—fl > 
1, ti + r* > 1. Then 

t^0(rt,ti) - ^^(r^fi) < {ti - fl)~^E^(ti) + {ti+rl)-^Ef{ti). 

Proof. This follows immediately from and the geometry of the region con¬ 
sidered. □ 

Proposition 11.11. With X as above, suppose 

ti<l.lto, |r*(ro)| + |r*(i7)| <.lto. (59) 


Then we have 


lf{X) < Ct-^E^ito). 


Proof. Let y be a smooth cut-off function equal to one on the interval [—1,1] 
and vanishing outside the interval [—1.5,1.5]. Define ^ to be a solution of the 
wave equation Ogtjj = 0 with initial conditions 

i’ito,r*) = x{‘2r*/to)(fito, r*), dt'il){to,r.,) = x(2r7to)9t<(>(to, r*). 

We claim that 

I^{X)=I^{X). 
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This follows from the choices (1^ and the fact that initial data for ' 0(^01 ’) 
coincides with those of 0(to, •) for the values of —.5to < r* < .Stg. 



We now apply Proposition to the function 0 to obtain 

4(7Z')<CE^(to). 


It remains to show that 


From the definition of the energy and the properties of the cut-off func¬ 
tion X we see immediately that 

.75to 

E^(to) < zz7^(-.75to,to) - w^(.75to,to) + C'io^ J j\(l)\‘^r^d(T dr*. 


-.75to §2 


In view of Proposition ll 1.1(11 it suffices to show that 
.75*0 


j j\ct>\\^dadr* <CE^{t^). 


-.75*0 §2 

We will rely on the one-dimensional inequality 


\f{x)\^dx <Ca^ ( j \dxf{x)\'^dx + j\f{x)\^dx 


Applying this to the function r0, and then integrating over S^, we obtain 
.75to / -TSto 

J J l^pr^dcrdr* < Ctl i J J r^dcrdr* 


.75*0 S2 


V.75*o S2 


i 

11 


-1 


< CEfito), 

where the last inequality follows from itCTl-ira. 


□ 
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12 Local observers’ uniform energy bounded¬ 
ness 

For any V 2 > 1 construct a rectangle TZ and triangle TZ as in the figure of 
Proposition im with vi = 1. Let 

CXD 

^<l>= J J r'^i^ - + {dr-(j)f + {I - {l,r* ,as2)dr*das2. 

— oo 

Similar to etc. we define etc. Note that < 

Eq. By Cauchy stability^^, we have that 

F^{[ui,oo] X {!}) < CH^. 

On the other hand, by Proposition IIl.i^ we have the uniform estimate 

I^{TZ) < CE^ 

and thus, by Proposition [TT31 we have 

(TZ) < C [E^ + H^), (60) 

F^{[ui,ao]x {v 2 }) < C {E^ +H^), (61) 

({oo} X [1, V 2 ]) < C{E^ + H^). 

Note: Using the above bounds, the classical result of Kay and Wald can be 
reproven without exploiting discrete isometries of Schwarzschild. See Section 

[TO 

13 Proof of Theorem 17.11 

Set to = 1, ti+i = l.lti, and let Ui, vt be defined by the relations r(ui,Vi) = ro, 
Ui + Vi = ti. Define the sets 

= [Mi,00] X [vi,Vi+i], 

K = {U<t< U+i}, 

s[ = yj)=o'Fl 

Xi =iz[f^ {ro <r < R}, 

y = 

see this, one must actually change coordinates. 
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13.1 Energy decay 

Applying first Proposition 111. HI to X, we obtain 

ifm < CE^, 

Ilm) < CE^^. 

Applying now Propositions 111.71 and 111.'^ with TV = S[, X = yi, we obtain 

E^[U)<CUE^,^^+CE^{1), 

and similarly 

Ef^[ti) < CUE^^,^_+CE^J1). 

Applying now Proposition lll.lll to Xi, we obtain 

< Ct-^E^^^^+CtfE^{l), 

and also 

iliX,) < + Ct-^El{l). 

Applying Proposition lll.il and once again Proposition lll.71 but now with TZ' = 
TZ'^, X = Xi, we obtain 

Summing over i, one obtains 

i 

j=o 

Proposition irn^ then gives 

E^iU) < C(logtdi^0.0„,0.„ + CEf^^Jl), 

and similarly, 

EliU) < C{\ogU)E^^,^^^,^^^^+CEl^^^Jl). 

One repeats the procedure one final time to remove the log term: 

EfiU) < (62) 

Note one also obtains then from Proposition lll.lll that 

If (Ad < CE,t-^. (63) 

Proposition III.KH (I62II and energy conservation now immediately imply m- 
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13.2 Local observers’ energy decay 

Refer to the following diagramme: 


O 



From 1(101) . we have already shown 

Apply Proposition 1 11. (II together with to obtain the existence of a Vi such 

that 

FJ ([iti, oo] X {?;}) < C{E^ + H^){vi+i - Vi)~^ < C{E^ + (64) 

Now let Ui be defined by Vi) = tq, and construct rectangle 7Zi and triangle 
7Zi, as depicted: 


O 



By Proposition nn and (EH), we have that 

w^{ui,Vi) - w^{oo,Vi) < FJ ([ui,oo] X {{ij) < C{E^ + E[^)t~^. 

On the other hand, by m and the fact that m ~ u ~ on 7?.^, we have 

- Wf|,{u^+l,Vi)) < CE2t~^- 
Finally, from (ESI, we have the bound 

ifiTZ.) < CE,t-\ 

Thus, we obtain from Proposition ITTKl and the inequality E^ + H^j, < Eq, the 
bound 

Fj{[ui,oo\ X {i^i+i}) < C{Ei 
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We now repeat the procedure to obtain 


Fj ([m„oo] X {uj) < C{E, + Eo)t-^. 


Further application of ProDOsition lll.l^ now yields 

([m,oo] X {?;}) < C{Ei +Eo)v~^, 


(65) 


F^{{u}x[v,v])<C{E^+Eo)v-\ 


( 66 ) 


Note that, by ProDOsition lll.4| this implies in particular, that if {u,v) is such 
that r(u,v) < vq and if {u,v) is such that r(u,v) < tq then 



\'^(j)\'^da§2dv < C{Ei + Eq)v ^ 


(67) 


'OO 


L 


das^du < C{Ei + Eo)v-‘^. 


( 68 ) 


Recall that these estimates are not available in this region from the usual 
energy estimate. 

More generally one easily shows the following 

Theorem 13.1. LetS be an achronal hypersurface in the closure of the exterior. 
Then FJ^(S) < C{Ei + Eo){v+{S))~'^, where v+{S) is as in the statement of 
Theorem o and FJ (S) denotes the flux ofTapY'^ through S. 

The above theorem applies in particular to subsets of the event horizon 

13.3 Uniform boundedness of 0 

We now reprove the classical result of Kay and Wald stated as estimate (ITTIl of 
Theorem rm First consider the region r > uq . We have 



By the Sobolev embedding on S^, 


OO 


OQ 



+ \dr»^'^(j){t, p,a;)P) r^(p) da§2dp 




Thus, for r > ro 


< Cr 


(69) 
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In the region r < rg we define u(v) so that r{u{v),v) = rg. Then 
\(t){u,v,uj)\ < \(j){u{v),v,uj) + j \du(j){u',v,uj)\du' 

J u{v) 

/ \ w , f f°° (.du4){u',v,uj)f ^ 

< \(j){u{v),v,uj)\+c[ / --- du 

\J U{v) ^ M y 

Since r{u{v)^ v) = ro, by (Pli we have 

On the other hand, by the Sobolev embedding, we have 


poo 
*/ u(v) 


[du(j){u' ,v,uj)Y 


1 - M 


du' <C 


’u{v) ^ d' 

+ {du^’^4’{u',v,u}))‘^) da§2du' 


{{du(l){u' ,v,uj)f + {du^(l){u' ,v,uj)y 


where the last inequality follows from Pronosition 11 1 .4l and (161 II . Thus, for 
r < rg, 

< C + H4>,4>u,,4'^u,) ■ 

Combining the estimates for r > rg and r < rg and observing that Eg = 

obtain the desired result. 

13.4 Pointwise decay for (f) 

13.4.1 |(/)| < decay near 

For a fixed large R consider the region 

U = {(u,u) S {r < .R} : \(j){u,v,u!)\ < v <v, r{u,v) < R}, 

for a C to be chosen later. Id is clearly open. Moreover, for sufficiently large C, 
Id O {r < R}r\{v <2} . We will show that Id is closed, and thus, coincides with 
r < R. For this, clearly it suffices to show 

\^\ < iC/2)v-^ (70) 


vuld {v > 2}. 

Let {u, v) £ldr\{v> 2}, and consider the null segment {u} x [fi— 1,0]. Note 
that this null segment is contained in r < i?. We have a Sobolev inequality 

V 

|<))p(u,u,w) < K J J + \dy(j)\'^+ \f(l>\‘^+ \fdycl)\'^ 

V-1 §2 

+ ITOP + I das2 dv, 
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in r < R. Our energy decay estimate and our local observers’ energy decay 
estimate in particular imply that, after commuting the equation twice 
with angular momentum operators O, 

V 

j J {\dv^\^ + + IWP + das 2 dv < E 2 i-\ 

V-l §2 


in r < i?. Choosing C sufficiently large, it remains to show, say that 

V / ^ \ 

( c‘^\ 

|(/)p(M, u, a;)(i(Ts 2 dv < —— 


8K 


V-l §2 


Let uo{v) be defined by r*(uo{v),v) = v — v + R* + l, in particular uo(u)+u = 

V — R* — 1 = t is independent of v. We have 

V V V2 U 

j J\cj)\'^{u,v,uj)da§2 dv = II \(j)\'^{uo,v,u!)da§2 dv + 2 J J J du 4 > 4 ’dcr§2du'dv 

v — 1 v — 1 ^0 

V V OQ ^ 

< I J \(j)\'^(uo,v)dag2 dv + J J J diT§2 du'dv 


V — 1 Uq 


V oo 

J J J {1 — da§ 2 du'dv. 


V — 1 Uq §2 

From the bounds EHi), we have 


|(/)f (uo(u),'!;)dcr§2 du < Ci^iu . 


v-l §2 


On the other hand, from our local observers’ energy estimate ra . we have 

CO 

{du(t)Y 


II 

uo §2 


1 - 


■das 2 du < C{Ei + Eq)v . 


Thus, choosing C sufficiently large, we are left with showing, say, the bound 

V CO / '' \ 

J J J {I-das 2 du'dv < v~‘^. (71) 


We have the inclusion 

[uo(u), oo) X [u — 1, u] C {(t, r*) : v — R* — l<t< oo, v — t—1 < r* < v — t}. 
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Thus, 


CX5 V — t 


J j j{1 — y)\(f)\^ das^du'dv < C J j j{1 — ^jL)\(j)\^d<7^2 dr* dt. 


We may write 

OO V — t 




v — R* — l-\-A v — t 


v — R* — l 


J J y (1 —/r)|(^p(icrs2 dr* dt = J J y (1 —/r)|(/)pdcrs2 dr* dt 

-1S2 5-t-lS2 

OO v — t 

+ I y y(l -/r)|(/)pdtTs2 dr* dt. 

The first integral can be estimated with the help of p|i . (EH), 

v — R* — l-\-A v — t 

j f f {I — fi)\(j)\^das 2 dr* dt < 24(1 + |r — ,R* — ID^^C 

< ACEi{R*fi-‘^. 


sup Ef {t) 


t^[v — R* — 1,?} —/?* — 1+A] 


(72) 


For the second integral we use the fact that, for A large enough, the region of 
integration above is contained in U, as r*{u, v) < R* — A + 1 and — 1 < r < {; 
there, so we have 

|</>| < Cv~^. 

On the other hand, we have (1 — ^) < he*' in this region. Thus 


OO V — t —A-\-R* -\-l 

,| 2 ^„ _ ^„* ^/V 2~-2 f dr* <hC'^e 


(1—/r)|(()pdcr §2 dr* dt < C'^v 


r* ^ ^-A+R*+l~-2^ 




-OO 


Choosing A to be sufficiently large, say A = 10 + logd + R* ^ and C large as 
before, and satisfying in addition, say 

C > 32AK{R*f {E2 + Ei+ Eq) , 

we obtain GD, as desired. Note that E 2 dominates Ei and Eq. 


13.4.2 Decay in r > i? 

We turn to r > R. First consider the region {r > d(} n {u > 1}. By the Sobolev 
inequality, 

r^\(^iu,v,u;)\'^ <C J |(^|Vdas2+C' J {\rf(^\'^ + Aff(j3\'^)r^das2. 

S‘2 S2 
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For any i? < f < i? + 1, for fc = 0,1, 2, and any r we have 

J\ri^(j)\‘^r'^{t,r,uj)da§2 < j \VL^(j)\^r^{t,f*,u})da§2 

§2 §2 

r* 

+ j J {\dr>n^(l)\\n’^(l)\+r-^\i^''(t>\^)r^dpdas2. 

f* §2 

Applying a pigeonhole argument in f and (isHi), m with we obtain that 
r can be chosen so as for 

J \n’^(l)\^r^{t,f*,u;)das2 < 

§2 

On the other hand, from (P|l. Pll. and (ESI applied to we have 


dpdCT§2 < i52rt 


rn S2 


r 

J J dpdas2 < E2rt~^u~^. 

rj S2 

Therefore, since in the region m > 1 we have that (say) t > r/2 , it follows that 
\(t>\ < E2{rtu)~^ < 2E2r~^u~^. 

On the other hand, since max(r, u) >1^/2 we also obtain 

\4>\<E2{i + \v\)-\ 

The region {r > ro} O {m < 1} O {t > 1} can be treated similarly by 
integrating out to spatial infinity using the estimates 

f f dpda§2 < E2 \u\~^, 

Jr* JS 2 

f [ \n^(j)\^r dpdcjs^ < E2r~^. 

Jr* Js‘^ 

In this region, one obtains 

\4>\ < Cr“^(l + \(j,\ < £’2?;“^(1 + |m|)“5. 

This completes the proof. 
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